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ABSTRACT 

The radio spectral index map of the Coma halo shows a progressive steepening 
of the spectral index with increasing radius. Such a steepening cannot be simply 
justified by models involving continuous injection of fresh particles in the Coma halo 
or by models involving diffusion of fresh electrons from the central regions. 

We propose a two phase model in which the relativistic electrons injected in the 
Coma cluster by some processes (starbursts, AGNs, shocks, turbulence) during a first 
phase in the past are systematically reaccelerated during a second phase for a relatively 
long time 1 Gyr) up to the present time. We show that for reacceleration time 
scales of ^ 0.1 Gyr this hypothesis can well account for the radio properties of Coma 
C. For the same range of parameters which explain Coma C we have calculated the 
expected fluxes from the inverse Compton scattering of the CMB photons finding that 
the hard X-ray tail discovered by BeppoSAX may be accounted for by the stronger 
reacceleration allowed by the model. 

The possibility of extending the main model assumptions and findings to the case 
of the other radio haloes is also discussed, the basic predictions being consistent with 
the observations. 

Key words: acceleration of particles - radiation mechanisms: non-thermal - galaxies: 
clusters: general - galaxies: clusters: individual: Coma - radio continuum: general - X- 
rays: general 



1 INTRODUCTION 

The intracluster medium (ICM) consists not only of the hot 
gas emitting thermal X-rays, but also of magnetic fields 
and relativistic particles. The existence of cluster magnetic 
fields and relativistic electrons is directly demonstrated by 
the presence of diffuse non-thermal radio emission detected 
in a number of clusters, the so called radio haloes. Further- 
more, Faraday rotation measurements toward discrete radio 
sources in clusters, coupled with the X-ray emission from 
the hot ICM, have indicated the presence of ~ fiG magnetic 
field strengths. 

The radiative lifetime of the relativistic electrons, 
mainly due to inverse Compton (IC) losses with the cosmic 
microwave background (CMB) photons, is relatively short 
(~ 10^^* yrs) so that the radio haloes suggest the presence 
of relatively fresh injected and/or reaccelerated electrons in 
the cluster volume (Jaffe 1977; De Young 1992; Tribble 1993; 
Brunetti et al. 1999). In general, models which assume a pri- 
mary origin of the relativistic electrons require continuous 



injection processes and predict synchrotron radio emission 
and IC emission in the EUV and hard X-rays (Ensslin et 
al. 1997, 1999; Sarazin 1999; Volk & Atoyan 1999), while 
secondary electron models skip the problem of the reacceler- 
ation (Dennison 1980; Blasi & Colafrancesco 1999) and pre- 
dict large gamma-ray fluxes from neutral pion decay which 
could be tested by future gamma-ray missions. 

Among halo sources. Coma C in the Coma cluster is 
the most famous and well studied example of diffuse radio 
emission in cluster of galaxies and can be used to test the 
different models of radio halo formation and evolution. 

At the lowest frequencies the radio emission of Coma C 
extends up to ~ 30 — 40 arcmin from the center (e.g. Cordey 
1985 at 151 MHz; Hanisch & Erickson 1980 at 43 MHz; 
Henning 1989 at 30.9 MHz); at higher frequencies sensitive 
radio images have been obtained at 327 MHz (Venturi et 
al. 1990) and at 1.38 GHz (Kim et al. 1990). By applying 
Gaussian fits it has been found that the 327 MHz FWHM 
(28x20 arcmin) is significantly larger than that at 1.38 GHz 
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(18. 7x 13.7 arcmin), but smaller than the low frequency size. 
The total radio spectrum of the Coma halo is steep (a ~ 
1.2) and the 327-1380 MHz spectral index map shows that 
the spectrum steepens rapidly with increasing distance from 
the center (Giovannini et al.l993). The steepening has been 
confirmed by 1.4 GHz data from the Effelsberg single-dish 
100 m telescope (Deiss et al. 1997). 

In principle, current models invoking a continuous in- 
jection of rolativistic electrons may explain the total syn- 
chrotron spectrum of Coma C, but they fail in reproducing 
the spectral steepening with radius (Ensslin et al. 1999; Volk 
& Atoyan 1999). Here it is important to realize that the dif- 
fusion velocity of the rolativistic particles is low in relation 
to their radiative lifetimes so that the spectral steepening 
cannot be related to the diffusion of fast ageing electrons 
from the central regions into the cluster volume (Berezinsky 
et al. 1997; Sarazin 1999). As a consequence, unless one is 
ready to accept the rather unphysical scenario of an injection 
modulated across the cluster volume, the spectral steepen- 
ing must be related to the intrinsic evolution of the local 
electron spectrum and to the radial profile of the cluster 
magnetic field strength. 

Additional information of great importance in con- 
straining the models of Coma C can also be obtained from 
observations at higher frequencies. Recently, Fusco-Femiano 
et al.(1999) have discovered an hard X-ray tail exceeding the 
thermal emission. Also in this case the existing models are 
not able to reconcile an IC origin of the hard X-ray flux with 
the estimates magnetic field intensities derived from RM 
measurements, the required field being at least one order 
of magnitude weaker than the RM estimates (e.g. Sarazin 
1999, Ensslin et al. 1999 and references therein). This has 
stimulated alternative proposals for the origin of the hard 
X-ray excess (Blasi, Olinto, Stebbins 2000; Ensslin 2000 and 
references therein). 

In view of these difficulties we investigate the alterna- 
tive possibility and potentiality of a two phase scenario for 
Coma C in which relativistic particles injected during a first 
phase arc then rcacceleratcd during a second phase possibly 
associated with a recent merger event £is indicated by optical 
and X-ray observations. A preliminary, much simplified ver- 
sion of the model has been presented elsewhere (Brunetti et 
al.l999). Here radiation and Coulomb losses and reacceler- 
ation efficiency arc fully taken into account as a function of 
the cluster radius and the assumption of a constant magnetic 
field throughout the cluster volume is released. A detailed 
description of the time evolution of the electron spectra is a 
distinctive feature of our model as compared to other current 
models (e.g. Sarazin 1999 and Ensslin 2000 for a review) 

The main aims of the present work are to: 

i) investigate the efl^ccts of energy losses and reaccel- 
eration gains on the energy distribution of the relativistic 
electrons and find out whether the radio spectral steepening 
may significantly constrain the space of parameters of the 
physical processes at work in Coma C; 

ii) obtain the radial trend of the magnetic field strength 
for the same range of the physical parameters in i) and com- 
pare it with the available information; 

iii) compare the derived IC emission of Coma C with 
the hard X-ray flux measured by BeppoSAX and find out 
whether the relatively high central values of the magnetic 
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Figure 1. The predicted observable (dashed lines) and intrinsic 
(solid lines) spectral index for the Coma halo are reported as a 
function of the distance from the center. The curves axe obtained 
from Deiss et al.(1997) relationships based on 327/1380 MHz 
(bottom dashed and solid lines at r=0) and 327/1400 MHz (single 
dish) observations. The points are taken from three slices of the 
spectral index map (327-1380 MHz) of Giovannini et al.(1993) 
across the halo center. 

field indicated by the RM results are compatible with the 
IC scenario. 

In Section 2 wc introduce the model; Section 3 is de- 
voted to a general presentation of the time evolution of the 
electron energy distribution subject to reacceleration and 
losses; the results concerning the radio properties of Coma 
C are given in Section 4, while those concerning the hard 
X-rays are given in Section 5. In Section 6 we discuss our 
findings and extend the model for Coma C to the case of the 
other haloes, while the main conclusions are summarized in 
Section 7. In the Appendix we present a more detailed dis- 
cussion of the theoretical aspects of the model. Hq = 50 km 
s~^ Mpc~^ is assumed throughout. 



2 THE TWO PHASE MODEL 

2.1 The radio spectral steepening 

By assuming a spherical symmetry of the relevant physi- 
cal quantities and a Gaussian spatial profile of the emis- 
sion coefficient, Deiss et al.(1997) have derived relationships 
between the observed and intrinsic spectral index distribu- 
tions of the synchrotron emission. In the case of the Coma 
halo they have shown that the FWHMs and flux densities 
of the radio observations (Venturi et al. 1990; Kim et al. 
1990) can be reproduced by a central intrinsic spectral in- 
dex a ^ 0.2 — 0.4, leading to a central observed spectral 
index a ~ 0.7, in agreement with Giovannini et al. (1993) 
findings. In Fig.l we plot the observed and intrinsic spectral 
indices predicted from Deiss et al. (1997) relationships based 
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on the measurements at 327 MHz (Venturi et al.l990), 1.38 
GHz (Kim et al. 1990) and 1.4 GHz (Deiss et al. 1997). Both 
the data at 1.38 GHz (obtained with a synthesis radio tele- 
scope) and those at 1.4 GHz (obtained with a single dish 
radio telescope), when combined with the 327 MHz data, 
show a marked steepening of the observed spectrum with 
distance from the center; this confirms that the steepening 
in the Giovannini et al.(1993) spectral index map is not due 
to flux density losses at 1.38 GHz but it is real. The points 
reported in Fig.l are taken from three slices across the halo 
center of the spectral index map (Giovannini et al. 1993). Al- 
though there is some scattering, the agreement between the 
observed and derived spectral distributions is fairly good. 

2.2 Injection Phase 

We investigate the application of a two phase model for 
Coma C. During the first phase fresh particles are continu- 
ously injected in the volume of the Coma cluster for a suf- 
ficiently long time. An injection phase in the past can be 
reasonably related to starburst activity, during which vio- 
lent galactic winds can accelerate relativistic particles (Zi- 
rakashvili et al. 1996; Ptuskin et al. 1997), and/or to the 
activity of AGNs producing a large amount of relativistic 
particles (Giovannini et al. 1993; Sarazin 1999). In addition, 
the high temperature of the ICM suggests that the ICM 
has passed through strong shocks, possibly during major 
mergers events when a significant portion of the shock en- 
ergy can be channeled into the acceleration of relativistic 
particles (Blandford & Eichler 1987; Jones & Ellison 1991). 
Moreover, it is also possible that during major mergers relic 
particles, released in the past from starbursts and AGN's ac- 
tivity, can be reaccelerated and injected at higher energies. 

The time evolution of the particle energy distribution 
during the injection, N{'y,t), is obtained by solving the ki- 
netic equation (e.g. Kardashev 1962): 



dN{-f,e,t) 

dt 



d_ 



N{j,e,t) 



+ g,„j(7,6',t) 



(1) 



with Qinj{'y,0,t) the injection function. 
Following Sarazin (1999; Fig.l) we assume that the cooling 
of the relativistic electrons is dominated by Coulomb (^) and 
radiation(7'^/3(z, 6))losses, the relativistic bremsstrahlung 
being negligible, so that 



dt ^ 



-/3(e,^)7 



(2) 



The Coulomb losses depend on the ICM density as (e.g. 
Sarazin 1999): 



1.2.10-^^n(l-f ^'^('^Z") 



75 



1.4- 10" 



(3) 



while the coefficient of the radiative losses is: 

l3{e,z)^1.9-W~^ {B'^sin^e + Bjciz)) s"^ (4) 

where 6 is the pitch angle, z the redshift and Bjc the equiv- 
alent magnetic field strength of the CMB. 
If the particles are isotropically injected, i.e. Qinj{'~f,d) = 
<5inj(7)/(47r), the only dependence of the electron spectrum 
on the pitch angle is given by the synchrotron loss term 
(Eq.^. If the pitch angle scattering is an efficient process 
then the electrons are continuously isotropized, the term 



B^sin^O should be replaced by 2/3B^ and the solution is 
isotropic, i.e. N{-f,e) = A''(7)/(47r). 

The IC losses depend on the redshift as {1 + z)'^ so that 
the cooling time of the relativistic electrons in clusters is ex- 
pected to be dominated by the IC process at moderately and 
high redshifts; this produces a rapid ageing of fresh particles 
injected at z > 0.3 (Sarazin 1999). Due to the presence of 
reacceleration that dominates the time evolution of the elec- 
trons during the second phase and since it starts less than 
~1 Gyr before the corresponding cluster age, in our model 
the dependence on z of the radiative losses is not important 
and has been neglected in the calculations. 
Since the injection phase is assumed to hold in the past for 
a time considerably longer than the the cooling time of the 
electrons, the energy distribution approaches a stationary 
spectrum (Sarazin 1999). Under stationary conditions and 



assuming a time-independent Qinj ~ Kef 
(|) yield: 



iV(7,e) 



Ke 



^ m 



Eqs.(m) and 



(5) 



which is (X 7 at low energies (where Coulomb losses 



dominate) and oc 7 
losses dominate. ^ 



^ at higher energies where radiative 



2.3 Reacceleration Phase 

When the injection of fresh particles is stopped or strongly 
reduced the energy distribution rapidly steepens and the 
number of high energy radio emitting electrons rapidly de- 
creases. 

In order to avoid the problem of the rapid ageing, we 
investigate the possibility that some reacceleration in the 
cluster volume is present during a second phase lasting less 
than about 1 Gyr up to the present cluster age. The presence 
of reacceleration mechanisms in cluster of galaxies has been 
pointed out by many authors studying the spectral ageing 
in head-tail radio sources (e.g. Jaffe & Perola 1973; Ekers 
et al. 1978; Fanti et al. 1981; Parma et al. 1999). 

As a simplification of the scenario, we assume that the 
second phase starts when the efficiency of the injection pro- 
cesses is strongly reduced or stopped. The general case in 
which the end of the injection and the start of the reaccel- 
eration do not coincide is given in Appendix A; the main 
results remain unchanged. The second phase may be related 
to a merger or a post-merger scenario in which shocks and 
cluster weather, powered by the energy released during a 
merger, can compensate at some level for the radiative and 
Coulomb losses. A number of authors have found evidence 
for recent merging in the Coma cluster (Briel et al. 1992; 
Vikhlinin et al. 1997). Although we discuss specifically the 
case of Coma C, the radio haloes are generally found in clus- 
ters with recent merger activity (e.g. Feretti & Giovannini 
1996); the general application of the two phase model to the 
radio haloes is discussed in Sect. 7. 2. 

If the observed radio spectral steepening is related to 
the intrinsic evolution of the local electron energy distri- 
bution rather than to an unphysical injection modulation, 



we will assume (5=2.5 throughout this paper 
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important constraints on the physical processes at work in 
Coma C can be immediately obtained. 

Indeed, the spectral steepening implies that the break 
of the synchrotron spectrum is in the range 0.3-1.4 GHz for 
a relevant fraction of the cluster volume. By making use of 
the basic synchrotron relationships this immediately yields 
a typical break energy of the electron spectrum at 75 ~ 
l'f{B^G)~^''^ and a radiative lifetime T ~ 5- 10^^ [Bf.aY^'^ s 
for relatively weak magnetic fields {B < Bjc ~ SfiG). 
Therefore, in order to generate Coma C the reacceleration 
during the second phase must essentially balance the ra- 
diation losses, i.e. the reacceleration efficiency should be 
X ~ l/T ~ 2 • 10"^^/i/B^G s"^ Moreover, a significant 
flattening of the relic electron distribution may only be at- 
tained for a reacceleration time r such that electrons of en- 
ergy 7 ~ 100 are accelerated to ~ 75 (Section 3 for details 
and Eq.(||)), so that r > ln{'yt/100)x~'^ ~ 0.6-0.7 Gyr. 

If the average magnetic field intensity depends on the 
radial distance r from the cluster centre, it clea rly follows 
that a reacceleration efficiency x('') oc l/-\/ B{r) would pro- 
duce a constant radio spectrum throughout the cluster. By 
allowing for a likely decrease of B with r, it then follows 
that a constant reacceleration efficiency results in a system- 
atic steepening of the synchrotron spectrum with r simply 
because at a given frequency higher energy electrons are 
selected by the lower field intensity. This steepening effect 
would be enhanced if the reacceleration efficiency increases 
in the regions of stronger fields. On the other hand. Coulomb 
losses are particularly severe in the central region of the clus- 
ter leading to a marked fiattening of the electron spectrum 
at low energies and, consequently, to a fiattening of the syn- 
chrotron spectrum at low frequencies as observed. 

The above qualitative description of the effects of sys- 
tematic reacceleration and losses on the evolution of the 
electron spectrum will be discussed in detail in the following 
Section 3. (Readers specifically interested in the application 
of the model to Coma C may skip this Section and go di- 
rectly to Section 4.) 



3 EVOLUTION OF THE ELECTRON AND 
SYNCHROTRON SPECTRA 

In this Section we illustrate the time evolution of the elec- 
tron energy distribution and related synchrotron spectra for 
a wide range of the relevant physical parameters. This pro- 
vides the basic quantitative scenario delimiting the region 
of the parameters' space which produces a self-consistent 
picture of the emission from Coma C. 



3.1 The electron spectrum 

As a simplification we assume that the particles are system- 
atically reaccelerated by first order Fermi mechanisms. The 
time evolution of the energy of the electrons is: 



J = -/3(e)72-^ + X7-'?(7)7 



(6) 



where x expresses the efficiency of the reacceleration and 
77(7) the relativistic bremsstrahlung losses given by (e.g. 
Sarazin 1999) 



X = 2.0*10 " s 




Figure 2. The predicted electron energy distribution is shown 
in arbitrary units as a function of the reacceleration period r 
as indicated in the panel. The calculations are performed from 
Eq.(|) for X = 2.0 ■ IQ-^'^ s'^ , ^ 



Given the typical cluster thermal density n ~ 
10~^cm— 3, the bremsstrahlung losses are always ~ 100 
times lower than the reacceleration efficiencies adopted in 
this work and will be neglected. 

With the initial {t = ti) spectrum given by Eq.(P) the 
solution of the kinetic equation (Eq.|l| with Q — Q) is: 



N{i,r,e) 



K,q{e) [l-tanh^(r^/2)]7-^+^ 
S-1 [276(r,e)tanh(r^/2)]2/33(e) 

S~3 / ^ _i \ -S+l 



7 



1-7/76(^,5) 



+ 



m J 



(8) 



where q{0) = x^-'i^f3{9) {q > 0; see Appendix A), r = t-U, 
and 76 is the break energy (i.e. the largest energy of the 
electrons) given by 



2m 



tanh^ 



(9) 



r;(7) ~ 1.5 ■ 10"^'^n [ln{-f) + 0.36] 



(7) 



For ^ = (no Coulomb losses) one finds 75(6') = xl PO- ~ 
e^^"^) in agreement with Kardashev (1962) results. 

An example of the evolution with time of the energy 
distribution of the electrons is represented in Fig. 2 for x = 
2.0 • 10~^^ s~^ (i.e. a reacceleration time ~ 1.6 x 10* yr) 
and 5 — 1.3 • 10~^^ s~^ (i.e. a thermal gas density n ~ 10~^ 
cm"^). 

Due to the radiative losses the high energy part of the 
distribution evolves rapidly with time. When the break en- 
ergy 7f, reaches the asymptotic value (x + \/— g)/2/3 the radi- 
ation losses are balanced by the reacceleration and the break 
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Figure 3. The predicted electron energy distribution is shown 
in arbitrary units as a function of the reacceleration efficiency as 
indicated in the panel The calculations are performed from Eq.(^ 
for T = 0.8 Gyr, ^ = 1.3- IQ-'^^s''^ and B = 0.0. 
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Figure 4. The predicted electron energy distribution is shown in 
arbitrary units as a function of the Coulomb losses as indicated 
in the panel. The calculations are performed from Eq.(^) for x = 
2.0 ■ 10-1'^ s-i, r = 0.8 Gyr and B = 0.0. 



energy is frozen; when r becomes sufficiently larger than the 
reacceleration time the combined effects of Coulomb losses 
and systematic reacceleration cause a depletion of the low 
energy side of the electron energy distribution which rapidly 
flattens. 

In Fig. 3 we report the electron energy distribution after 
a representative time interval r = 0.8 Gyr for different reac- 
celeration efficiencies. As expected, one finds that stronger 
reaccelerations yield higher energy breaks and flatter energy 
distributions. 

The effect on the particle energy distribution due to 
the increase of the Coulomb losses is illustrated in Fig.4 for 
a given reacceleration. By increasing the fCM density the 
Coulomb break in the initial spectrum is positioned at larger 
7's, causing a larger flattening of the electron spectrum at 
lower energies (Eqj^, so that the reaccelerated spectrum is 
also flatter. At this point it should be noticed that, since 
the ICM density falls down with increasing distance from 
the center, also the electron energy distribution is expected 
to be flatter in the central regions of the clusters. 

The effect of increasing the radiative losses is illustrated 
in Fig. 5 for a given reacceleration. As expected, the balance 
between losses and reacceleration is obtained at lower ener- 
gies with increasing B. Furthermore, with increasing B the 
break energy 'yb becomes closer to the energy at which the 
reacceleration balances the Coulomb losses and the energy 
distribution slightly flattens. 

All these findings hold in general: in the case of sys- 
tematic reacceleration of relativistic electrons diffuse in a 
thermal medium (as in clusters of galaxies) after a time suf- 
ficiently large, which depends on the reacceleration and/or 
losses, the electron energy distribution is expected to be flat 
or slightly inverted before a sharp cut-off. 



3.2 The synchrotron spectrum 

In general the synchrotron emissivity is obtained by inte- 
grating the synchrotron Kernel over the electron energy and 
angular distributions (cf. Westfold f959; Pacholczyk f970). 
We assume that the electron momenta are isotropically dis- 
tributed and that the fields are tangled and of constant in- 
tensity on a sufficiently small scale. These assumptions allow 
to consider the synchrotron emissivity (averaged over a suf- 
ficiently large volume) isotropic. By omitting for simplicity 
of notation the dependences of the functions, the emissivity 
per unit solid angle is : 



^ Ut' V 5-1 4 mc2 



r/2 



desine4>{e) 



v?sin9 



duF\ Jd:±-A{i-uf-''u-'+^ 



1 + 



2,l + (Cw-^-X7i>(e))/(7?(W))i2 , i 



-5+1 



1 - u 

where the electron Kernel 

|>oo 

F{y) = y Kdv)dv 



+ 



(10) 



K& being the 2/3 order modified Bessel function. 



u{e) = 



1 -tanh2(r^/2) 
[tanh(r^/2 ; 



-(4+2)^-3 



p-%0) 



= c 



B'^sin^e + Bjc 



Ihie.r) ^+ y^/tanh(l^) 



(11) 

(12) 
(13) 
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Figure 5. The predicted electron energy distribution is shown in 
arbitrary units as a function of the magnetic field strength B as 
indicated in the panel The calculations are performed from Eq.(^ 
for X = 2.0 ■ 10-1*^ s-l, ,e = 1.3 . IQ-I^ s"!, 9 = 55° and r = 0.8 
Gyr. 
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Figure 6. The predicted synchrotron spectrum is shown in ar- 
bitrary units for two well separated values of the Coulomb losses 
as indicated in the panel. The frequency scale is in unit of the 
break frequency i/j,. We assume two reacceleration efficiencies: 
X = 2.5 • 10-1*5 s-l (solid lines) and X = 4.2 • lO'^'^ s"! (dashed 
lines). The calculations are performed for a reacceleration period 
T=0.8 Gyr and B = l.O/xG (typical for the Coma center). 



C ~ 3 • 10 ^ cgs and i^j, the critical frequency calculated for 
7 = 76 (0 = 90°). 

As it is well known, the shape of the synchrotron spec- 
trum depends on the electron distribution in the 3-D mo- 
mentum space; this is a function of the combined losses 
and reaccelerations (Eq.^). The IC losses are constant over 
the cluster volume, but both the magnetic field B and the 
Coulomb losses are expected to vary with the distance r 
from the center. Thus, the relativistic electron energy dis- 
tribution (Eq.^, the break energy and the importance of the 
synchrotron losses with respect to the IC losses depend on 
r. As a consequence the shape of the synchrotron spectrum 
(Eqs.|l^, |l^ depends on r and it is expected to be inter- 
mediate between those calculated in the case of pure syn- 
chrotron losses (Kardashev 1962; Pacholczyk 1970; Komis- 
sarov & Gubanov 1994 : KP-model) and of pure IC losses 
(Jaffe & Perola 1973; Komissarov & Gubanov 1994: JP- 
model) . 

As discussed in the preceeding Section increasing the 
Coulomb losses and/or the systematic reacceleration yields, 
in general, a flattening in the energy distribution of the rel- 
ativistic electrons. The effects on the synchrotron spectrum 
are illustrated in Fig. 6 for given r and B. As expected, the 
synchrotron spectrum below the break frequency becomes 
flatter with increasing the ICM density and/or the system- 
atic reacceleration. It should also be noticed that the syn- 
chrotron specra are not sensitive to the presence of spectral 
spikes in the electron energy distribution (Fig.3). If stochas- 
tic Fermi mechanisms provide a substantial contribution to 
the reacceleration of the electrons, the spectral spikes are 
expected to be smoothed out, but our results remain un- 
changed. 



4 MODEL RESULTS FOR COMA C 

As shown in the previous Section a natural consequence 
of the combined reacceleration and Compton losses is the 
flattening of the electron spectrum and of the related syn- 
chrotron emission in the central region of a cluster, in qual- 
itative agreement with what is actually found in Coma C. 
In this Section the physical parameters will be constrained 
in order to reproduce the observed spectral steepening and 
total radio spectrum. 



4.1 Basic assumptions 

Following Deiss et al.(1997) we approximate the brightness 
distribution at each frequency with a symmetric Gaussian 
of FWHM A(t/) equal to the geometric mean of the smaller 
and larger FWHM measured at each frequency. The very 
rapid decrease of the radio brightness given by the Gaussian 
approximation cannot well constrain the synchrotron elec- 
tron distribution at large distances from the cluster center 
where the brightness fades into the noise. Thus, in the model 
computation, the tail of this distribution should be allowed 
to vary within the observational uncertainties. To this end 
we model the radio brightness distribution by adding to the 
Gaussian from the fits at each observed frequency a symmet- 
ric Gaussian of larger FWHM A;(i/) so that the observed 
emissivity is represented by: 
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Aiy) giving the normalization of the radio brightness dis- 
tribution and /, a free parameter (from the radio maps we 
find / < 0.06). 

As shown in the preceding Section the Coulomb losses play 
an important role in shaping the electron energy distribution 
and the associated synchrotron spectrum. By adopting the 
/3-model of Briel et al.(1992) the Coulomb losses are : 



^(r) =4.2- 10 



-3/3/2 



(15) 



where J?c=10.5 arcmin is the core radius of the ICM and 
(3 = 0.75 (the central density is no ~ 3 • 10~^ cm~"^). 

We assume that the efficiency of the reacceleration can 
be parameterized as the sum of a constant large scale com- 
ponent and a small scale component. Indeed, turbulent gas 
motion originating from the massive galaxies orbiting in the 
cluster can amplify the efficiency of the large scale reacceler- 
ation within the cluster core (Deiss & Just 1996) . As a pos- 
sible choice, we assume that the small scale reacceleration is 
proportional to the inverse of the typical distance between 
galaxies. By adopting the observed galaxy distribution in 
the Coma cluster (Girardi et al. 1995), we parameterize the 
reacceleration as : 



X{r) = XLS + XG. 



(16) 



Rg=5.5 arcmin being the optical core radius and aa^O.S. 

It should be noticed, however, that since the radio 
brightness and the spectral properties of the model are ob- 
tained by integrating the emissivities over the cluster vol- 
ume, the main model results do not strongly depend on the 
precise choice of the functional dependence in Eg. (p^ . 

In the sequel we adopt a reacceleration time interval 
r = O.SGyr, but we anticipate that the essential features of 
the model are not modified as long as r is of the order of 1 
Gyr. 

At last we assume complete reisotropization of the elec- 
tron momenta, not only for simplicity, but also because the 
pitch angle scattering is likely to be an efficient process dur- 
ing the reacceleration phase. We have carried out several 
tests and found that, due to the low magnetic field intensi- 
ties, the differences between the synchrotron spectrum (and 
all the other relevant quantities in this paper) calculated in 
the isotropic and anisotropic case are negligible. 



Matching the observations: the inferred B 
field distribution 



4.2 



From Eqs. (|1C|-|16D , at each r and for a given set of the param- 
eters (xls, Xgo, and B), we obtain the synchrotron emis- 
sivity as a function of the frequency measured in terms of 
the break frequency i'i,{r) (vbir) oc B{r)j^{r; 9 — 90°) ) and 
then we can compute the spectral index between 327 MHz 
and 1.4 GHz. 

In Fig. 7 we plot the calculated intrinsic spectral index 
between 327 MHz and 1.4 GHz at r = 5 arcmin as a function 
of B and for different values of the reacceleration efficiency 
X- 

Small values of B lead to very steep radio spectra because 
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Figure 7. The calculated 0.327-1.4 GHz spectral index is given 
as a function of the magnetic field _B at a distance r = 5 arcmin. 
The calculations are performed for several reacceleration param- 
eters X, as indicated in the panel, and for the value of ^ obtained 
from Eq.(|l^). The reacceleration time interval is r = 0.8 Gyr. 



the break frequency of the synchrotron spectrum is compa- 
rable or smaller than 327 MHz. Moderate values of the mag- 
netic field allow to obtain very flat spectra as those derived 
for the central regions of the Coma halo (a\\ ~ 0.4 - 0.5 
at r = 5'). Since from simple synchrotron relationships the 
break frequency is Vi, cx B{B\q -I-B^)"^, it is found that the 
value of B for which the spectral index is the smallest for any 
given reacceleration x is obtained when B ~ Bic — 3^G, 
while the value of B for which one obtains a given spectral 
index depends on x- 

The magnetic fleld B{r) is numerically derived by im- 
posing that the calculated spectral index matches the in- 
trinsic one of Fig.l (the average between Deiss et al (1997) 
and Giovannini et al (1993) curves). The results are shown 
in Fig. 8 for relevant values of the reacceleration parameters. 
For relatively small distances from the cluster center the 
value of the magnetic fleld is sensitive to the assumption 
about both xls and xg, while at larger distances, due to 
the rapid decrease of the small scale reacceleration efficiency 
(Eqjl5[), the value of B depends only on Xls- 

With reference to Fig. 7 and a = 0.5, we find a lower 
bound to the central value of the reacceleration efficiency 
x(0) ~ 3.3 • 10~^® s~^ and a corresponding magnetic field 
strength B{0) ~ 3ij.G which is the largest value of the central 
field intensity compatible with our model. 

We emphasize that these results depend only on the 
radio spectral indices and not on the normalizations of the 
relativistic electron distributions. 

The magnetic field strength required by our model to 
match the spectral indices in the central region can be com- 
pared with the values given in the literature, obtained by 
Faraday rotation measurements. Kim et al. (1990) have es- 
timated a magnetic field B = (1.7 ± 0.9)/^G' in the central 
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Figure 8. The B distribution required by the model to match 
the spectral observations is reported for different reaccelerations 
as indicated in the panel. The calculations are performed for a 
reacceleration period t = 0.8 Gyr. Due to the reacceleration pro- 
file (Eq.p^) the magnetic field in the majority of the cluster vol- 
ume depends only on XLS- The dots represents the Jaffe (1980) 
theoretical predictions for the magnetic field distribution in the 
Coma cluster (we have assumed Bq = 3fiG). 



regions of the Coma cluster, assuming a tangling scale of 
the magnetic field of 10-40 kpc. More recently Feretti et 
al. (1995), studying the RM distribution of the head-tail 
radio source NGC 4869 located near the center of Coma, 
inferred the existence of a strong magnetic field component 
{B ~ 6fj,G) tangled on the kpc scale. The two values are 
not in conflict, considering that the 6 fiG field could be as- 
sociated to the local environment around NGC4869. The 
value obtained by our model in the central Coma region ( 
~ 1 — 3/xG) is consistent with the above measurements. It 
is also in agreement with the results of Clarke et al.(1999), 
who, using statistical RM technique in a sample of low red- 
shift rich Abell clusters, find central values of the magnetic 
field of the order of a few fiG. 

The RM results can be used to constrain the reacceler- 
ation parameter from above. With reference to Fig. 7 we find 
that the central value of the magnetic field is significantly 
smaller than 1 fj,G for x(0) > 5- 10~^^ s~^, setting an upper 
bound to the reacceleration efficiency in our model. 

The origin and structure of the magnetic fields in galaxy 
clusters are presently debated. Seed fields could be generated 
by turbulent galactic wakes (Jaffe 1980) and/or by ejecta 
from galaxies (Kronberg et al. 1999; Volk & Atoyan 1999) 
and maintained by turbulent dynamo driven by galaxy mo- 
tions through the ICM (De Young 1992). Such low level mag- 
netic fields can also be amplified by cluster mergers (Nor- 
man & Bryan 1998; Eilek 1999). According to Jaffe (1980), 
the magnetic field distribution depends on the thermal gas 
density (n) and on the distribution of massive galaxies 
(nc) as B(r-) = B(0)(n(r)/n(0))"-'(nG(r)/nG(0))°-'*. Jaffe 
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Figure 9. The total synchrotron radio spectrum predicted by 
our model is compared with the flux measurements for the three 
model's parameters specified in Fig. 8: model a) (dotted lines), 
model b) (dashed lines) and model c) (solid lines). The calculations 
are given for two different values of / (Eq.^): 0.01 (lower lines), 
0.02 (upper lines). 

(1980) predictions calculated for the Coma cluster, adopting 
B{0) ~ S^iG and with n(r) from Brief et al.(1992) and nair) 
from Girardi et al.(1995), are consistent with our model find- 
ings (Fig.8). 

4.3 The predicted total synchrotron spectrum 

The model at this point is anchored to the observed and in- 
trinsic 327-1400 MHz spectral index distribution as a func- 
tion of r and to the observed FWHMs of the brightness 
distribution at these frequencies. 

The spectral shape of the synchrotron emission at a 
given radius is very different from a simple relatively steep 
power law such as the total Coma C radio spectrum. The 
shape of the electron energy distribution changes with r in 
a predictable way depending on losses and gains which, as 
shown in Section 4.2, are now well constrained. 

As a consequence, a fundamental test of the model is 
whether it can reproduce the observed total synchrotron 
spectrum. This is obtained by integrating the synchrotron 
emissivity (Eq.lO) over the cluster volume taking into ac- 
count the r-dependence of the parameters involved in the 
calculation. It is at this point that the electron spectra are 
normalized as a function of r in order to match the ob- 
served emissivities (Eq.l4) derived from the brightness dis- 
tributions at 327 and 1400 MHz following Deiss et al (1997) 
procedure. As an illustrative example in Fig. 10 we plot the 
electron energy distribution at different distances from the 
cluster center for the representative case c) of Fig.8. 

The calculated total synchrotron spectra for the same 
set of parameters of Fig.8 are represented in Fig.9 and com- 
pared with the observations. Once the uncertainties in the 
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Figure 10. The electron energy spectra are reported at different 
distances from the cluster center (normalization scale in arbitrary 
units). The calculations are performed for model c) with /=0.015 
(Fig. 8). The energy interval of the electrons emitting hard X-rays 
via IC scattering of the CMB photons is also shown (HEX). 

observations are considered, the model well accounts for the 
observed radio spectrum over all the frequency range sam- 
pled by the observations. Since the model is mainly forced 
on the 327 and 1400 MHz observations, it is not surprising 
that the computed spectra between these two frequencies are 
essentially indistinct. On the contrary, at higher and lower 
frequencies the calculated synchrotron spectrum is more sen- 
sitive to the assumed parameters. In order to see whether 
this effect can help in further constraining the model pa- 
rameters it is important to assess the quality of the spectral 
measurements. 

The radio fluxes quoted in the literature are obtained 
by subtracting from the observed fluxes those contributed 
by point sources falling within the halo. Since the obser- 
vations at different frequencies have been performed with 
very different instruments and, in some cases, at times dif- 
fering of many years, the subtraction of the sources is not 
standard and can influence the calculation of the halo flux. 
In particular, the low frequency points at 43 and 74 MHz 
in Fig. 9 were obtained by subtracting the extrapolated flux 
from 5C4.81 and 5C4.85 and a contributed flux from the 
unresolved sources measured at 430 MHz and extrapolated 
down to low frequencies (Hanisch & Erickson 1980); the 
30.9 MHz halo flux has been obtained by subtracting all 
the sources from the 151 MHz survey with an extrapolated 
flux density at 30.9 MHz larger than 1 Jy (Hanning 1989). 
More precise subtraction procedures have been applied at 
higher frequencies thanks to the better spatial resolution: 
the 327 MHz radio flux of the Coma halo has been obtained 
by subtracting the 327 MHz flux densities from a list of 64 
unresolved sources observed at 327 and/or 1400 MHz (Ven- 
turi et al. 1990) and a similar procedure has been adopted 
by Deiss et al.(1997) at 1.4 GHz by using the list of ra- 



dio sources as compiled by Kim et al.(1994) (the 2.7 GHz 
flux is likely to be underestimated (Deiss et al. 1997) ). We 
find that the detailed high frequency subtractions are deeper 
than those performed at lower frequencies by 20-40 %, so 
that the 30.9, 43, and 74 MHz fluxes reported in Fig.9 may 
be considered as upper limits. It follows that the low fre- 
quency data cannot discriminate among the the models of 
Fig.9. 

In our model the radio spectrum at low frequencies 
< 100 MHz) is dominated by the contribution of the 
radiation emitted at large radii (r > 12 — 17 arcmin). The 
contribution from the particle with very low break frequen- 
cies at r > 30 — 35 arcmin, not detected in the high fre- 
quency maps due to ageing, is also important. Larger values 
of /, corresponding to an increased number of particles at 
large distances from the center, do not appreciably modify 
the 327-1400 MHz spectrum, but increase the low frequency 
radio fluxes expected by the model (Fig.9). 

By increasing the large scale reacceleration (xls), the 
spectrum of the synchrotron electrons is depressed and the 
expected radio fluxes at low frequencies decrease (Fig.9). On 
the other hand, by decreasing xls the low frequency radio 
spectrum steepens and may significantly overproduce the 
observations (Fig.9). This allows us to obtain a range for the 
large scale reacceleration efficiency of 1.5- 10~^^s~^ < Xls < 
3-10~^^s~^. Accurate measurements at low frequencies, such 
as those obtainable with the VLA at 74 MHz, are of crucial 
importance to constrain the model parameters. 

Due to the very steep radio spectrum of the external 
regions, the high frequency spectrum of the Coma halo does 
not depend on / (Fig.9). It is emitted by regions at r < 12 
arcmin and depends on both the large scale and small scale 
reacceleration efficiencies. With increasing reacceleration in 
the cluster core the particle spectrum fiattens and the break 
energy increases: this causes a decrease of the required B 
value so that the 327 MHz emitted frequency is closer to 
the synchrotron break frequency. As a net consequence the 
high frequency radio spectrum steepens (Fig.9). 

4.4 Relic electrons and halo energetics 

The shape of the injection rate (projected) of the relic elec- 
tron population, Ks{r), is shown in Fig. 11. It has been ob- 
tained by inverting the synchrotron emission equations for 
model (c) of Fig. 8 (the integral of Eq^| along the line of 
sight). The shape is rather close, within a factor ~ 2, to 
both the thermal and galaxy number density distributions 
in Coma. This result is very encouraging since no quan- 
titative assumptions on the space distribution of the relic 
electrons have been made in the derivation of the model. 
However, the similarity with the matter distribution is con- 
sistent with our conjecture that relativistic particles must 
have been efficiently injected into the ICM during the first 
phase by AGN and galaxy activity. 

The moderate upturn of Ke at r > 45 arcmin is mainly 
caused by the flattening in the assumed Coma radio bright- 
ness due to the superposition of the two Gaussian functions 
(Eqjl4[) . In any case this upturn is not relevant for the main 
characteristics and predictions of the model. 

The total energy associated with the relativistic elec- 
trons in the cluster is the result of the energy injected dur- 
ing the first phase and of the energy gained by the electrons 
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Figure 11. The projected distribution of the injection rate of 
the relativistic electrons Ke{r) (Eq.H) during the first phase is 
given (soUd lines) as a function of r (in arcmin) for model (c) of 
Fig. 8 (/= 0.01 and 0.02 starting from the bottom). The shapes 
of the projected density of the thermal gas (obtained from Briel 
et al.l992 fit) and of the galaxies (Girardi et al.l995) (dashed 
lines) are conveniently scaled for a closer comparison; the pro- 
jected galaxy density profile is only valid up to 38 arcmin from the 
cluster center. The projected radial distribution of the UV galax- 
ies (Donas et al.l995), normalized to the Girardi et al. counts, is 
also plotted for r > 25 arcmin. The region in which the 80-90 % 
of the IC hard X— ray flux is produced during the second phase 
(Section 5) is enclosed by the vertical dotted lines. 



during the second phase. We find that for a reacceleration 
parameter x ~ 3 ■ 10~^'^s~^ the total energy in relativistic 
electrons is ~ 10^^ erg. 

We point out that this value is much smaller (by a factor 
> 10) than that associated with the relativistic electron pop- 
ulation considered in the other proposed models (Kempner 
& Sarazin 1999 and references therein). This is because the 
rather flat electron energy spectrum, characteristic of our 
model, is much more efficient than those of the other mod- 
els in emitting synchrotron radiation at radio wavelengths. 



5 THE HARD X-RAY TAIL IN THE COMA 
CLUSTER 

From a deep BeppoSAX observation Fusco-Femiano et 
al.(1999) discovered an hard X-ray tail, in the range 20- 
80 keV, exceeding the extrapolation of the thermal X-ray 
emission spectrum; this excess has also been confirmed by 
an RXTE observation (Rephaeli et al. 1999) with consistent 
results. Due to the poor statistics, the spectrum is not well 
constrained so that the origin of this high energy emission 
cannot be firmly established. Several models have been ex- 
plored. It can be generated by the IC scattering of the radio 
emitting electrons with the CMB photons (Fusco-Femiano 
et al. 1999; Sarazin et al. 1999; Sarazin & Lieu 1998); in this 



case models invoking a secondary production of the elec- 
trons yield a gamma-ray fiux considerably larger than the 
EGRET upper limit (Blasi & Colafrancesco 1999). Alter- 
natively, it might originate via relativistic bremsstrahlung 
of a supra-thermal power law tail of particles (Ensslin et 
al. 1999), or, as proposed by Dogiel (2000), emitting from 
the thermal emission of a modified maxwellian distribution 
of the hot ICM induced by strong acceleration processes 

The IC origin is particularly attractive in view of the 
fact that a comparison between the radio emission from the 
halo and the hard X-ray tail discovered by BeppoSAX, if 
spatially coincident, may allow a direct estimate of the av- 
erage strength of the magnetic field B over the cluster vol- 
ume. Nevertheless, since in the present model the electron 
energy spectrum depends on the distance from the cluster 
center (Fig. 10), the synchrotron as well as the IC spectra 
are emitted by relativistic electrons whose energy distribu- 
tion cannot be represented by a unique power law over the 
whole volume. Therefore, the magnetic field cannot be sim- 
ply evaluated by standard formulae (e.g. Harris & Grindlay 
1979). 

Under our assumptions (Eq.^) the IC emissivity per 
unit energy and solid angle in the ultra-relativistic Thomp- 
son approximation is given by: 
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e being the energy of the CM B pho tons and for ultra- 
relativistic electrons = 61 /4e (e.g. Blumenthal & 
Gould 1970). As a cluster volume we assume a sphere of 
radius 55 arcmin (~2.2 Mpc), that is about the radius up 
to which the counts of UV galaxies, presumed cluster mem- 
bers, begin merging the background galaxies (Donas et al. 
1995). 

Since the IC emission is only sensitive to the number 
of electrons, it is expected to be mainly contributed by the 
outer parts r > 30 of the cluster volume which contain the 
large majority of the relativistic electrons able to IC scatter 
the CMB photons in the BeppoSAX hard band (Fig. 10). As 
a consequence, the predicted X-ray luminosity is sensitive to 
the large scale reacceleration efficiency; we find that the from 
80 to 90% of the model 20-80 keV flux is emitted between 
30-50 arcmin from the center (Fig. 11). 

The calculated hard X-ray IC fluxes from our model 
for the three choices of the parameters given in Fig.8 are 
compared with the BeppoSAX data in Fig. 12. Models with 
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Figure 12. The IC fluxes from the Coma cluster for the three 
models of Fig.8 (lower line /=0.01, upper line /=0.02). The Bep- 
poSAX hard X-ray data are obtained from Pusco-Femiano et al. 
(1999) with the spectral indices at the 90% confidence level; a 15% 
uncertainty in the normalization is assumed. The cross represents 
the Coma EUV excess Bowyer et al.(1999). 

relatively high reacceleration efficiencies may account for a 
large fraction or all of the X-ray flux measured by Bep- 
poSAX. By comparing Fig. 12 with Fig. 9 it is seen that it is 
precisely this type of models that contribute less to the syn- 
chrotron flux at very low radio frequencies. It is not surpris- 
ing, then, that in order not the overproduce the BeppoSAX 
flux we obtain an upper limit for the large scale reaccelera- 
tion efficiency xls < 3 ■ 10~^®s~^, the same as that found 
in Section 4.6. Obviously, accurate measurements at low ra- 
dio frequencies would be of crucial importance not only to 
constrain the parameters of our model for Coma C, but also 
to check whether an IC origin of BeppoSAX hard X-rays is 
viable. Thus, at variance with statements found in the lit- 
erature, central magnetic fields of ~ l/j-G may be consistent 
with an IC origin of the BeppoSAX hard X-ray flux, the 
reason being that while the X-rays are mainly produced in 
the outer regions of the cluster the radio flux above a few 
100 MHz is mainly produced in the inner regions. 

In Fig. 12 we have also represented the EUV flux from 
Bowyer et al.(1999): the model EUV flux is well below the 
observed. The origin of the EUV excess is still unclear, its 
spatial distribution is considerably narrower than the radio 
brightness distribution and, if non-thermal, it is probably 
associated to an additional electron population (Bowyer & 
Berghofer 1998; Ensshn et al. 1999; Brunetti et al. 1999). 



6 DISCUSSION AND CONCLUSIONS 

6.1 A two phase model for Coma C 

It has been shown that the radio properties of Coma C can 
be explained by a two phase model consisting of an injection 



period, in which a power law spectrum of rclativistic elec- 
trons has been continuously injected throughout the cluster 
volume and modified by radiative and Coulomb losses, fol- 
lowed by a reacceleration pheise of the aged (relic) spectrum 
for typical duration of ~1 Gyr up to the present time. 

The model is anchored to the spectral index map of 
Coma C between 327 and 1400 MHz which shows a marked 
steepening going in the outward direction from a flat central 
value. The explanation of this pattern, not accounted for in 
other published models, constitutes the basic motivation of 
this work. 

On the basis of simple considerations (Section 2.3) we 

have shown that the reacceleration efficiency should be x 
2 • lO^^^s^^ and the reacceleration time r > 0.6 Gyr. A 
detailed modelling has then been worked out based on the 
following assumptions: 

• The brightness distributions at 327 and 1400 MHz 
are described as the sum of two gaussians; the radio images 
constrain the peak amplitude of the large scale gaussian to 
be a fraction / < 0.06 of that of the small scale gaussian. By 
assuming a spherical symmetry, the emissivities at 327 and 
1400 MHz and the corresponding spectral index as a func- 
tion of the radial distance r arc obtained from the observed 
brightness distributions by making use of Deiss et al. (1997) 
procedure. 

• The Coulomb losses, ^(r), are flxed by the /?-model 
flt of the ICM of Briel et al (1992). 

• The reacceleration efficiency is assumed to be rep- 
resented by the sum of two terms: a large scale reaccel- 
eration (xls) constant throughout the cluster volume and 
a small scale reacceleration (xG{r)) peaking at the cluster 
center. We have identified the large scale component with 
shocks and/or turbulence, possibly generated during a re- 
cent merger, and the small scale component with the am- 
plification of such shocks/turbulence by the motion of the 
massive galaxies in the cluster core. 

• Let evolve the relic electron population, subject to 
radiation and Coulomb losses, for a reacceleration time r = 
0.8 Gyr. 

With these assumptions we find that: 

• The electron energy distribution depends on r, be- 
ing much fiatter and with a somewhat higher energy break 
toward the cluster center where both the magnetic field in- 
tensity and the reacceleration parameter are stronger and 
the Coulomb losses are more severe. 

• The value of the 327 1400 MHz spectral index in the 
central region (a = 0.5) implies an upper bound to the in- 
tensity of the magnetic field B ~ 3/uG and a lower bound to 
the reacceleration efficiency x(0) ~ 3.3 ■ 10^^^' s^^. 

• A central magnetic field of at least 1 )iG, as indicated 
by RM measurements, implies an upper bound to the reac- 
celeration efficiency at the center of x(0) ~ 5 ■ 10~^®s~^. 

• By fitting the steepening of the 327-1400 MHz spectral 
index with r, the magnetic field intensity smoothly tapers 
off from central values of 1-3 fiG to ^ 0.1/iG at r ~ 50 
arcmin. 

• The model synchrotron emission perfectly fits the to- 
tal radio spectrum of Coma C within 327 and 1400 MHz, in- 
dependent of the X and / parameters allowed to vary within 
the limits set above. This demonstrates the internal con- 
sistency of the model. On the contrary, the computed low 
frequency fluxes vary appreciably with the parameters, be- 
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ing higher with decreasing x and increasing /: most of this 
flux is emitted by the outer parts of the cluster volume, so 
that by a closer comparison with the observed low frequency 
flux it is possible to set limits on the large scale acceleration 
efliciency 1.5 ■ 10^^'' < xlS < 3 • IQ-^^s'^. 

• With the electron energy and space distributions con- 
strained by the radio data it is found that the X-rays pro- 
duced by IC scattering the CMB photons can account for 
the hard X-ray flux discovered by BeppoSAX if the larger 
flgures of the xlS and / parameters, allowed by the radio 
model, apply. Like the flux emitted at low radio frequen- 
cies, most of the X-rays are produced in the outer volume of 
the cluster. This demonstrates once again the crucial impor- 
tance of accurate flux measurements at frequencies < 100 
MHz. We remark that the BeppoSAX flux is signiflcantly 
overproduced for xls > 3 ■ 10~^®s~^ (/ = 0.02), which is 
consistent with the upper limit set by the radio spectrum at 
low frequencies. 

• By de-evolving the electron spectra constrained by the 
radio model it is found that the shape of projected radial 
distribution of the injection rate of the relic electron pop- 
ulation is very similar to that of the galaxies in the Coma 
cluster. This is consistent with our conjecture that the relic 
electron population may have been injected by AGN and 
galaxy activity at earlier epochs. 

• The present energy in relativistic electrons required 
by our model is ~ 10^^ ergs, at least a factor 10 smaller 
than in other proposed models. The energy budget required 
at the end of the injection phase is ~ 2 • 10^* ergs, if al- 
most immediately followed by the reacceleration phase, and 
becomes larger by increasing the time gap between the two 
phases simply because the spectrum of the relic electrons 
is further depressed. We find that a temporal gap of 3-4 
Gyr involves an energy budget ~ 100 times larger, but still 
reasonable (see Appendix). 

We emphasize that the results listed above arc essen- 
tially obtained by imposing that the model accounts for 
the radio properties of Coma C and we find that this re- 
quirement is already stringent enough to constrain the rele- 
vant parameters, such eis the reacceleration's efficiencies and 
time, within a narrow range. Wo also notice that the typical 
reacceleration efficiency, corresponding to an e-folding time 
scale of ~ 0.1 Gyr, is relatively modest and comparable or 
smaller than the values frequently quoted both in clusters 
(cluster weathering) and in tailed radio galaxies. 

Finally, the likely presence of relativistic protons in- 
jected with the electrons deserves some attention. It can 
be constrained by looking for the gamma-ray emission from 
the decay of the tt" generated in the p-p collisions. Recently, 
Blasi & Colafrancesco (1999) have worked out a secondary 
electron model for Coma C and estimated an upper bound 
- 4 ■ lO^^erg {Ho = 50 km s"^ Mpc"^) for the total energy 
of centrally injected protons in order not to exceed the up- 
per limit on the gamma ray flux sot by EGRET (this figure 
should be much larger in our model where the majority of 
the relativistic particles are injected in regions of low ICM 
density). Since in our model the particle reacceleration is 
due to shocks and/or turbulence in the ICM, the energy 
gained by relativistic protons during the second phase can 
not exceed the thermal energy of the ICM (~ lO^'^erg) so 
that the gamma-ray flux from p-p collisions in our model 
is expected to be far below the EGRET limit. 



6.2 General considerations for the radio haloes 

Although the present work is focussed on Coma C, we be- 
lieve that a two phase modeZ may explain the radio properties 
of other haloes as well. 

Like for Coma the required energy input in the ICM 
could be supplied by a recent merger as suggested by the ev- 
idence of a significant association between radio haloes and 
cluster mergers. Moreover, it has been found that the occur- 
rence of radio haloes in clusters increases with the cluster 
X-ray luminosity Lx in the 0.1-2.4 keV band (from ~ 4% 
for Lx < 5 ■ lO^^erg s"^ to ~ 36% for Lx > 10''"''erg s"^ 
Giovannini 1999 in preparation). This is consistent with the 
expectations of the two phase model since more X-ray lumi- 
nous, massive clusters would have produced a larger number 
of injected particles during the first phase, thus supplying 
sufficient energy for the formation of the radio haloes also 
after a moderately long gap between the injection and the 
reacceleration phases. 

It should be born in mind, however, that the strength 
of the halo radio emission will critically depend on the com- 
bination of the various parameters, (reacceleration efficiency 
and duration, magnetic field strength and the temporal gap 
between the two phases) so that we expect that only a frac- 
tion of the massive clusters with a recent merger will show a 
substantial radio halo. For instance, contrary to expectation 
a too eSicient reacceleration may result in an electron spec- 
trum very much stretched toward higher energies so that, 
other conditions being equal, there may be fewer relativistic 
electrons emitting at radio frequencies in a given magnetic 
field and the halo emission would be dimmed. 

The IC emission of the relativistic electrons with the 
CMB photons will typically fall in the UV and EUV band at 
the end of the injection phase, progressively moving toward 
the X-ray band during the reacceleration phase. Thus, non- 
thermal X-ray emission from clusters with recent mergers 
should be a quite common feature, depending on the total 
number of relic electrons. In general, as seen in the case of 
Coma, the radio at intermediate frequencies and the hard 
X-rays arc expected to be mainly contributed by difi^crcnt 
regions : the former near the center of the cluster, associated 
with regions of relatively higher B, the latter from the outer 
parts of the cluster volume. Clearly, observations with future 
hard X-ray facilities will be of great value in checking the 
correctness of this scenario. 
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APPENDIX A: THE GENERAL TWO PHASE 
MODEL 

In Section 3 we have analyzed the effect due to the sys- 
tematic reacceleration on the energy distribution. In order 
to simplify the scenario, it was assumed that the reacceler- 
ation phase starts when the injection phase is stopped or 
strongly reduced. 

In this Appendix we generalize the results of this paper 
by assuming that the two phases are separated in time. 

After the injection is stopped at a given time ti, we let 
the particles to age for a time interval At. Mostly due to 
IC losses with the CMB photons, the particles rapidly age 
and the energy distribution is depopulated at relatively high 
energies. 

By integrating Eq.(j^ with constant losses one obtains: 



7(e,Af) 



1 + 7(t.)\/m7Ctan(v^^At) 



and a break energy 

lb Ad) 



Ifili tan(ypAt) 

Then, from the kinetic equation (Eq|l|; with Q 
finds the energy distribution of the electrons : 



(Al) 



(A2) 



0), one 



iV(7,At,e) 



Ji'e[l + tan2(^At)] 



/3(e)(5-l) 



7 

76,1 



7 + 



-ihb. 



(A3) 



which becomes Eq.(|^ for At 0. 

We assume that after a time At following the end of 
the injection phase, the particles are reaccelerated. The time 
evolution of the energy of the electrons is obtained by inte- 
grating Eq.(pl): 



g < ^ 
g > 



X-2e7(^) 
79 



tan 



-1 x-2S7(At) 



tanh-i x^2Ml) 



tanh' 



-1 x-2e7(At) 



(A4) 



where q{9) ~ — 4,^fi{9) expresses the importance of the 
reacceleration. 



Eq.(A4) can be explicated in terms of 7; it gives 
^ 2^7(At)-x- V~9tan(a;) x 



^^^'^^ 2^ (257(At)-x)tan(x)-v 
for g < 0, while 

7(Af)(yg/tanh(x)+x)-2^ 



i{r,e) 

for g > 0, 

l^l 



2/3(e)7(At) + ^/tanh(a;)-x 



(A6) 



(A7) 



and ''" is the reacceleration period. Eq.(A.6) reduces to 
Eq.CT) for r = 0. 



The break energies are obtained from Eqs.(A.5 & A6) by 
imposing 7 (At) 00 : 




3 

Log 7 



Figure Al. The predicted electron energy distribution, in the 
same arbitrary units as in Fig. 3, is shown as a function of the 
reacceleration period r as in dicated in the panel. The calculations 



are performed from Eq. ( A9 ) for x ■ 
and ; 



: 2.0-10" 



-1,^ = 1.310- 



: 0.23. The parameters A an d 7t, 1 can be calculated 
respectively from Eqs.{All) and (A15); one has : 3.4 Gyr, 230 



(solid Une), 3.24 Gyr, 238 (dotted linoy72.9 Gyr, 255 (short dashed 
line) and 2.6 Gyr, 271 (large dashed line). 



g < ^ 

^{2etan(a;)}-^ + x/2C 
g > ^ 

{X + V9/tanh(x)}(2/3(e))-i 



(A8) 



In the case of the Coma cluster, g > 0; this is because by 
assuming a physically sound reacceleration time ~ lO^yr it 



is X 



10" 



while, given the thermal gas density and the 



IC losses, it is 4^/3 < 4 x 10~ . The energy distribution of 
the reaccelerated electrons is obtained by so lvin g the kinetic 
equation (Eq.|l} with Q = 0) with Eqs.(^, ^ & 



A(7,r,e) = 



q{e)Ke l-tanh-2(x) 1 + t&n^ {^/ l3{e)iAt) 



s-i 



(■5+1) 



2 c 

m 



im - 



where 



m = 1 + 



ib{9)m 



m-ib 



ib{e) 



{i-ihmr (A9) 



(AlO) 



The final shape of the energy distribution of the emit- 
ting electrons is mostly dependent on the reacceleration 
phase during which it is heavily modified. For this reason, 
even if the particles can be injected and age over cosmolog- 
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ical time scales, in the calculation of Eq.(A3) we have not 



considered the variation with z of the IC losses. 

However, the final energy distribution, Eq.(A9), de- 
pends on At and on 7i,_i which are directly related to the 
redshift of the end of the injection phase and to the redshift 
at which the reacceleration phase starts; such dependences 
can be simply introduced in the model. 

In the case of an Einstein-de Sitter universe, one has 
that : 



(All) 



Zi is the redshift corresponding to the end of the injection 
phase, while Zacc, the redshift corresponding to the start of 
the reacceleration phase, is given by: 

-2/3 



SHoT 



(A12) 



If the particles age over cosmological time scales, the IC 
losses are expected to dominate the Coulomb losses for 7 > 
100 so that Eq.(^) can be simplified and solved by consider- 
ing the dependence of /? on 2 and that dz/dt = Ho\/l + z ; 
one finds: 



l{Zac 



where 



b{zi, 



1 + ^b{z^ 



(A13) 



(A14) 



the break energy is simply obtained from Eqs.(A13- A14) 
with 7(«i) 00: 



76,1 



x5/2 



{l + Za 



-.5/2 - 



(A15) 



Given r and Zi, from Eqs .( |A12| [A15| , All) one has 
76,1 and At, then from Eq.(A9) one has the particle energy 



distribution at 2 = 0. 

In Fig.Al the model energy distribution of relativistic 
electrons injected up to a redshift 2 = 0.23 is shown for 
different reacceleration times. The break of the electron en- 
ergy distribution rapidly increases for a reacceleration time 
r < 0.8 Gyr, then it reaches the asymptotic value ~ 10*. 
During the reacceleration the IC emission from such a pop- 
ulation is expected to be channeled at the beginning in the 
UVs, then in the soft and in the hard X-rays. 

By assuming that the energy released during a merger 
generates an efficient reacceleration of the order of ~ 10*yr 
in the cluster volume for ~1 Gyr, then non-thermal hard 
X-ray emission should be produced for ~ 0.5 Gyr. 

The possibility of a gap between the two model phases 
gives a strong limit in the detection of non-thermal emis- 
sion (either synchrotron and IC) from merging clusters. De- 
pending on the gap duration, a considerable fraction of the 
electrons is thermalized or shifted at low energies (where 
the Coulomb losses dominate also during the reacceleration 
phase) after the injection phase. 

We find that, if the gap is ~ 3 — 4 Gyr long, the energy 
content of the relativistic electrons is reduced by a factor ~ 
100 (Figs.3 and Al). 
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